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\ The energy of a graph is defined as the sum of the absolute values of the eigenvalues of the graph. In this 



o 



X 



Abstract 



paper, we present a new method to compare the energies of two /c-subdivision bipartite graphs on some cut 
edges. As the applications of this new method, we determine the first L^^j^J largest energy trees of order n 
as : for n > 31, and we also give a simplified proof of the conjecture on the fourth maximal energy tree. 
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1. Introduction 



Let G be a graph with n vertices and A be its adjacency matrix. Let Ai, • • • , A„ be the eigenvalues of A, 

n n n 

then the energy of G, denoted by E(G), is defined [2, N] as E(G) = J^M. 



i=i 



> 

The characteristic polynomial det(a;7 — A) of the adjacency matrix A of a graph G is also called the 
characteristic polynomial of G, written as <fi(G,x) = a*i(G)x n ~ % . 

i=0 

In this paper, we write bi(G) — \ai(G)\, and also write 



4(G,x) = ][>(G) a 



If G is a bipartite graph, then it is well known that </>(G, x) has the form 



LfJ LfJ 
<t>(G, x)=J2 ^(G)x n - 2t = J2(-iyb 2l (G)x n - 21 (1.1) 



and thus 

LfJ 



0(G, x)=Y, b 2l (G)x n ~ 2 \ (b 2t (G) = \a 2l (G)\ = (-l) 4 a 24 ) (1.2) 



i=0 



In case G is a forest, then b 2 i(G) = m(G, i), the number of i-matchings of G. 
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The following integral formula by Gutman and Polansky (Q) on the difference of the energies of two 
graphs is the starting point of this paper. 



E(Gi)-E(G 2 ) 



hi 



cj>(Gi,ix) 



Ax 



(i 



(1.3) 



4>(G 2 ,ix) 

— oo 

Now suppose again that G is a bipartite graph of order re. Then by and (ll.2[) we have 

<t>{G, ix) = i n ${G, x) (G is bipartite, i = ^/-l) (1.4) 

Using (|1.4I) we can derive the following new formula from (| 1 . 3[) which does not involve the complex 
number i. 

Theorem 1.1. If G\,G 2 are both bipartite graphs of order re, then we have 



E(G X )-E(G 2 



mt^d, 

<KG 2 ,x) 



(1.5) 



Proof. Since G\, G2 are both bipartite graphs of order n, it is easy to see that 

LfJ 



0(Gi,x) _ j=o 



£ 6 2j -(Gi). 



r n ~ 2 j 



MG 2 ,x) LfJ 
So from (|1.3[) and (|1.4|) we have 



is an even function and ^L- — — — - > for x > 0. 

0(G 2 ,x) 



C(G X )-E(G 2 ) =- 



hi 



cj)(G u ix) 



<j>{G 2 ,ix) 



+00 



dx = 



In 



4>{G x ,x) 



d>(G 2 ,x) 



+ OO 



dx = 



, ^{Gux) A 

m — da:. 

4>(G 2 ,x) 



□ 



Definition 1.1. Let f(x) — Yl a-iX n ~ % and g{x) — biX n ~ l be two monic polynomials of degree n with 
nonncgative coefficients. 

(1) . If a. L < hi for all < i < n, then we write f(x) =4 g(x). 

(2) . If f(x) =4 g(x) and f(x) ^ g(x), then we write f(x) -< g{x). 

Now we define the following quasi-order for bipartite graphs (which is equivalent to the well known 
quasi-order defined by the coefficients bi(G) ). 

Definition 1.2. Let Gi and G 2 be two bipartite graphs of order n. Then we write Gi =^ G 2 if (f>(Gi,x) =4 
4>{G 2 ,x), write Gi -< G 2 if (f>(Gi,x) -< (f>(G 2 ,x) and write Gi ~ G 2 if (f>(Gi,x) — <j){G 2 ,x). 

According to the integral formula in Theorem ll.il we can see that for two bipartite graphs Gi and G 2 of 
order re, 

Gi ^G 2 =>E(Gi) <E(G 2 ); and Gi -< G 2 =>■ E(Gi) < E(G 2 ). 



The method of the quasi-order relation "=^" is an important tool in the study of graph energy. 



2 



Graphs with extremal energies are extensively studied in literature. Gutman [if determined the first and 
second maximal energy trees of order n; N.Li, S.Li Q determined the third maximal energy tree; Gutman et 
al. Q conjectured that the fourth maximal energy tree is P n (2, 6, n — 9) (see Figj3]for this graph); B. Huo 
ct al. 7] proved that this conjecture is true. 

In this paper, we first consider in §2 some recurrence relation of the polynomials cf>(G(k), x) for the k- 
subdivision graph G(fc) (on some cut edge e of a bipartite graph G). Then in §3 we present a new method of 
directly comparing the energies of two fc-subdivision bipartite graphs G{k) and H(k) if they are quasi-order 
incomparable. Using this new method, we are able to provide a simplified proof of the above mentioned 
conjecture on the fourth maximal energy tree. The main result of this paper is that, we determine (in §5) 
the first L^ir^J l ar § es t ener gy trees of order n > 31 by using the new method of comparing energies given in 
§3. For example when n > 2007, we can determine the first 1000 largest energy trees of order n (but up to 
now, only the first four are known). 

2. Some recurrence relations of 4>(G, x) and <fi(G, x) for fc-subdivision bipartite graphs 

The following lemma is an alternative form of Heilbronner's recurrence formula 0. 
Lemma 2.1. jjg/ Let uv be a cut edge of a graph G, then 4>(G, x) = <p(G — uv, x) — <j>(G — u — v, x). 




For the sake of simplicity, we sometime abbreviate <j>(G,x) by <fi(G). 
The following relation can be derived from Lemma 12.11 

Lemma 2.2. Let Hi, H2, H3 be graphs as shown in Fig\]\ Then we have 

<j)(H 3 , x) = x(f>{H 2 ,x) — <j>{Hi,x) 

Proof. Let G[ be the graph obtained from G\ by attaching a new pendent edge uw to G\ at u, and G' 2 be 
the graph obtained from Gi by attaching a new pendent edge vy to G 2 at v. Then by using Lemma 12.11 we 
have 

<l)(G' 1 )=x<t>(Gi)-(l>{G 1 -u), and 4>{G' 2 ) = x4>(G 2 ) - <j>(G 2 - v). 
Now using Lemma 1 2. II for H3 and its cut edge e' = wy, we have 

<t>{Hi) = 0(^3 ~ e) - <f>(H 3 -w-y) = HG'i)HG 2 ) - 0(Gi)0(G 2 ) 
={x<t>{Gi) - <P(Gi - u))(x<f>(G 2 ) - <P{G 2 - v)) - 0(Gfi)0(Ga) 

= (x 2 - 1)0(G!)0(G 2 ) - xc/>(G 1 )cf > (G 2 -v)- ^(G 2 )0(G 1 - u) + <f>(Gi - u)<j>(G 2 - v) 
Also using Lemma \2. II for H 2 and H 2 — e\ we have 

cb{H 2 ) =0(i? 2 - ei) - (f)(H 2 ~u~w) = <f>{H 2 -e x - e 2 ) - <j>{H 2 -e 1 -w-v)- 0((Gi - u) U G 2 ) 
=xcj>{G 1 )cj>{G 2 ) - 0(G!)0(G 2 - v) - - u)<f>{G 2 ) 
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Using Lemma |2~T1 for Hi we also have 

tf>{H{) = <t>{Hi - e) - tf,{Hi -u-v) = 0(Gi)^(G 2 ) - 0(d - u)0(G 2 - v) 
Now it is easy to verify from the above three equations that (j)(H 3 ) = x0(77 2 ) — <j)(H\). 



□ 




G G{k) 

Fig. 2: Graph G and its fc-subdivision graph 



Definition 2.1. Let e be a cut edge of a graph G, and G e (k) denote the graph obtained by replacing e with 
a path of length k + 1 (for simplicity of notations, we usually we abbreviate G e (k) by G(/c) ). We say that 
G(k) is a fc-subdivision graph of G on the cut edge e. We also agree that G(0) = G. 

From Lemma [2~2l we have the following recurrence relation for cf>(G(k),x). 

Theorem 2.1. Let G(k) be a k- subdivision graph of G on the cut edge e of G, then we have 

4>{G(k + 2), x) = xcf>(G(k + l),x) - <p{G(k),x) (k > 0) 

Proof. Take Hi = G(k) in Lemma 12.21 and e be an edge in Hi on the path of length k + 1 obtained by 
fc-subdividing the edge e. Then H2 = G(k + 1) and H3 = G(k + 2). The result now follows from Lemma 
l2~2l □ 



Theorem 2.2. Let G be a bipartite graph of order n and G(k) be a k-subdivision graph (of order n + k) of 
G on some cut edge e. Then we have 



(j>{G{k + 2), x) = x<p{G(k + 1), x) + <p(G(k),x) (k > 0) 

Proof. By Theorem 12.11 we have 

(j){G{k + 2),x) = x<t>{G{k + 1), x) - 4>(G(k), x) 

substitute x by ix, we get 

4>(G(k + 2), ix) = ix(j){G{k + 1), ix) ~ <f>(G(k),ix). 
Now using (|1.4[) for G(k + 2), G(k + 1) and G(k) (since they are all bipartite) we have 
i n+k+2 4>{G{k + 2), x) = i n+k+1 ix0(G(k + l),x) - i n+k 0(G(k),x) 
Dividing both sides by i n + k + 2 we get (|2.ip . 



(2.1) 



□ 
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Theorem 2.3. Let e, e' be cut edges of bipartite graphs G and H of order n , respectively. If G(0) =5! H (0) 
and G(l) =^ H(l), then we have G(fc) =4 H(k) for all k > 2, with G(fc) ~ H(k) if and only if both the two 
relations H(0) ~ G(0) and H(l) ~ G(l) hold. 

Proof. The result follows directly from Theorem 12.21 and induction on k. □ 

Theorem 2.4. Let G,H be bipartite graphs of order n, ei,e% be two cut edges of G and e^e^ be two cut 
edges of H . Let G(a,b) denote the graph obtained from G by subdividing e\, e.<i by a,b times, and H(c,d) 
denote the graph obtained from H by subdividing e^e^ by c,d times, respectively. If 

G(0,0H #(0,0) and G(0, 1) ==! H(0, 1), (2.2) 

(7(1,0)=$ #(1,0) and G(l, 1) =$ H(l, 1) (2.3) 

t/ien we have G{l,k) =4 H(l,k) for all I > and k > 0. Moreover, if one of I and k is at least 2, then 
G(l,k) -< H(l,k) if each of \2.2\) and h2.3\) contains at least one strict relation. 

Proof. Using Theorem 12.31 for ei and e' 2 we have 

HO I =>> G(0, fc) ^ ff (0, k) (k > 0), (2.4) 
(1231) G(l, fe) ^ i?(l, fe) (fc > 0). (2.5) 

Now using Theorem 12.31 for ei and e[ we also have 

(2.4) and (2.5) => fc) ^ fl"(i, jfc) (I > 0). 
When (]2.2[) and (]2.3[) both contain strict relations, we have both strict relations in (|2.4|) and f|2 . 5[) for 
k > 2. Thus G(Z, fc) -< i?(/, fc) for all k > 2 by Theorem 12.31 Similar arguments apply to the case I > 2. □ 



3. A new method of directly comparing the energies of fc-subdivision bipartite graphs 

Notice that if the conditions in Theorem 12.31 do not hold, then G(k) and H(k) might be quasi-order 
incomparable. In this section, we present a new method to directly compare the energies of two fc-subdivision 
bipartite graphs G(fc) and H{k) when they are quasi-order incomparable. Using this method, we give a 
simplified proof of the conjecture on the fourth maximal energy tree. 

In the following, we always write gk = 4>(G(k),x), hk = 4>(H(k),x), and dk = ——■ 

9k 

Lemma 3.1. Let G{k), H{k) be k-subdivision graphs on some cut edges of the bipartite graphs G and H of 
order n, respectively (k > 0), gk, hk and dk be defined as above. Then for each fixed x > 0, we have 

(1) . If d\ > do, then do < dk < d\ for all k>2; 

(2) . If di < d , then d% < d k < d for all k > 2; 

(3) . If di = do, then dk = do for all k. 
(So in any case we have dk > min{d , di}.) 

Proof. By the recurrence relations in Theorem 12. 2\ we have 

h k xhk-i+hk-2 xd k -i9k-i + dk-29k-2 

d k =— = : = : 

gk xg k -i + gk-2 xg k ~i + gk-2 

= ( XJ ¥ ) d ^ + ( gk ~2 ) d *-> 

\xgk-i+gk-2j \xg k -i+gk-2j 
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This tells us that dk is a convex combination of dk-i and dk-2 with positive coefficients, which implies 
that dt lies in the open interval (dk-i, dfc-2) or (dfc-2, dfc-i) if dk-i ^ d/c-2- Using this fact and the induction 
on fc we obtain that dk always lies in the open interval (do, d\) or (di, do) when do ^ d\, and dk = do when 
d\ = do- n 

The following theorem can be derived from Lemma |3. II 

Theorem 3.1. (1). If hxg - h^gi = 4>{H(l),x)4>(G(0),x) - (f>(H(0), x)(f>(G(l), x) > (which is equivalent 
to di(x) > d{j{x)) for all x > 0, then we have 

E(H(k)) - E(G(fc)) > E(iT(0)) - E(G(0)) (for all k > O.J 

(2). Ifhig -h g 1 = <p(H (1), x)<p(G(Q), x) — <p(H(Q), x)(j>{G{l), x) <0(whichis equivalent to d x (x) < d (x)) 
for all x > 0, then we have 

E(F(fc)) - E(G(fc)) > E(if(l)) - E(G(1)) for allk^l. 

Proof. (1). Since di(x) > do(x) for all x > 0, by (1) of Lemma \3. II we have dk(x) > do(x) for all x > and 
fc > 0. So by (fL~5|) we have 

+00 



E(F(fc))-E(G(fc)) = i / ln gg(fe)^) dx = £ /" i n4(a; ) da; 

TT J 6<G(k),x) I" J 



4>(G(k),x) 


+00 +00 — 

>- / lnd (a;)da; = - / In tf g ^' ^ da: = E(H(0)) - E(G(0)) (fc > 0). 
n { n { 4>(G(0),x) 

The proof of (2) is similar to that of (1). □ 

In [9], Shan et al. show that the fourth largest energy tree is either P„(2, 6, n — 9) or T n (2, 2|2, 2) (see Fig|3] 
and Fig|4]for the definitions of these two graphs). B. Huo et al.[7] proved that the conjecture on the fourth 
maximal energy tree is true by showing that E(P n (2,6, n — 9)) > E(T n (2, 2|2, 2)). Now by using Theorem 
13.11 we are able to give a simplified proof of the conjecture on the fourth maximal energy tree. 

Theorem 3.2. If n > 10, then 

E(P„(2,6,n-9))>E(T n (2,2|2,2)) 

Proof. Let H = Plo(2, 6, 1) and G = Tio(2, 2|2, 2), e be the pendent edge on the pendent path of length 1 in 
H, and e' be the edge between the two vertices of degree 3 in G. Then we have P n (2, 6, n — 9) = H(n — 10) 
and T„(2, 2|2, 2) = G(n — 10). By some directly calculations, we have 

(f(H(0), x) = <MPi (2, 6, 1), x) = x w + 9 x 8 + 27 x 6 + 31 x 4 + 12 x 2 + 1, 
0(G(O), x) = 4>{T l0 {2, 2|2, 2), x) = x w + 9 a; 8 + 26 x 6 + 30 x 4 + 13 x 2 + 1, 
4>{H{l),x) = 4>{P U {2, 6, 2),x) = x 11 + 10 a; 9 + 35 x 7 + 52 a; 5 + 32 x 3 + 6 x, 
4>{G{l),x) = 4>{T lx {2,2\2,2),x) = x 11 + 10 a; 9 + 34 x 7 + 48 a; 5 + 29 x 3 + 6 x. 

So we have 

<j)(H(l),x)<i>(G(0), x)-$(H(Q),x)$(G(l), x) = 2a; 15 +22a; 13 +89a; 11 +168a; 9 +156a; 7 +66a; 5 +9a; 3 > (x > 0). 
Also by using computer we can obtain 

E(H(0)) = 11.937511, E(G(0)) = 11.924777, So E(H(0)) - E(G(0)) = 0.012734 > 0. 
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So by Theorem l3.il we have for n > 10, 

E(P n (2, 6, n - 9)) - E(T„(2, 2|2, 2)) = E(iJ(n - 10)) - E(G(n - 10)) > E(H(0)) - E(G(0)) > 0. □ 

Combining Theorem 13.21 with the result that the fourth largest energy tree is either P„(2,6, n — 9) or 
T„(2, 2 1 2, 2) ([9J]), we conclude that the fourth maximal energy tree is P n (2, 6, n — 9). 

Remark: Here we would like to mention that, the main points of the simplification in the proof of Theorem 
GUare: 

1. We use the integral formula (|1.5|) (instead of f|l .3[) ) which uses the real polynomial <fi(Gj,x) instead of 
the complex polynomial 4>(Gj, ix) for j = 1, 2. 

2. The recurrence relation (|2.1I) for <fi(G(k),x) allows us to use Lemma I3TT1 to directly compare dk{x) and 
do(x) (namely directly compare the integrands \ndk(x) and m<io(x) in the formula (11.51) for K(H (k))— E(G(fc)) 
and E(H(0))— E(G(0))), without the need of solving the recurrence relation (|2.1[) to obtain explicit expressions 
for h k = <fr(H(k), x) and g k — (j)(G(k), x). □ 

Notice that in Theorem l3.1[ we need either di(x) > do(x) for all x > or do(x) > di(x) for all x > 0. Now 
if both of these two conditions do not hold, then both d (x) and di(x) are not a lower bound for c4(a;) (k > 2). 
Although in this case we can not use Theorem 13. 1[ but by Lemma 13.11 we still have min{d (a;), di(x)} as a 
lower bound for dk(x) (for all x > 0). Thus we can still have the following lower bound (which is independent 
of k) for E(iJ(fc)) -E(G(fc)). 

Theorem 3.3. Let G(k), H(k) be k-subdivision graphs of bipartite graphs G and H on some cut edges. Let 
4>(H(k), x) 

dk{x) = — ■ — and let D — {x > 0|<io(a;) > ^i(^)}. Let D c be the complement of D in (0,oo). Then : 

4>{G{k),x) 

+ 00 

MH(k)) - MGVk)) > - j \nm\n{d {x),d l {x)}dx = - [ In d x (x) dx + - [ \nd {x)dx (3.1) 

IT J 7T J 7T J 

o D D c 

where the right hand side of A3.1\) can also be written as: 

+oo 



— [ \ndi (x) dx H — f In do (x) dx = — [ \ndi (x) dx f In d\ (x) dx H — f In do (x) dx 

7T J 7T J IT J 7T J IT J 

D DC DC D c 

=E(H(l)) E(G(1)) - - / In ^f\dx (3.2) 

7T J d (X) 
DC 

or equivalently, 

- [ lnd 1 (x)dx + - [ \ndo(x)dx = M(H(0))-R(G(0)) + - [ In^ldx (3.3) 

IT J IT J TT J do(X) 

D D c D 

Theorem 13.31 will be used several times in §4 and §5 in the proof of our main results. 

4. Some upper bounds for the energies of non-starlike trees 

In the following discussions, we will divide the trees into two classes. One is called the starlikc trees, and 
the other one is the non-starlike trees. In this section, We will give some upper bounds for the energies of 
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the non-starlike trees. We will show that the energy of a non-starlike tree is bounded above either by the 
energy of P n (l, 2,n- 4), or by the energy of T n (2, 2|2, 2) (see FigEJand Figgjl. 

Let Ns(G) be the number of vertices in G with degree at least 3, and A(G) be the maximal degree of G. 
A tree T is called starlike if N 3 (T) < 1, and is called non-starlike if N 3 (T) > 2. 

It is easy to see that if N 3 (T) = 0, then T is the path P n . Now if N 3 (T) = 1, then T consists of 
some internally disjoint pendent paths starting from its unique vertex with degree at least 3. Suppose 
that the lengths of these pendent paths are positive integers a\,a,2, • • • , <2fe. Then we denote this tree T by 
P n (ai, a-2, ■ ■ ■ , o,k), where a± + a.2 + ■ ■ ■ + = n — 1 and k = A(T) (see Figj3j). Sometimes we also denote 
P n (ax, a.2, ■ ■ ■ , Qfe) by P„(ax, a?, • • • , ctfc-i, *), since * is uniquely determined by n and oi, 02, • • • , flfc-i- 




Let a, b, c, d be positive integers with a + b + c + d < n — 2. Let T„(a, 6|c, d) be the tree of order n obtained 
by attaching two pendent paths of lengths a and b to one end vertex of the path P n -a-b-c-d, and attaching 
two pendent paths of lengths c and d to another end vertex of the path P n -a-b-c-d (see FigHJ. 

It is not difficult to see that if T is a tree of order n with A(T) = 3 and N 3 (T) = 2, then T must be of 
the form T n (a, b\c, d), where a + b + c + d<n — 2. 




Fig. 5: The graph G u (a,b) 



In [9| and [10(, Shan et al. studied how graph energies change under edge grafting operations on unicyclic 
or bipartite graphs and proved the following result in the comparison of the quasi-order on unicyclic or 
bipartite graphs: 

Lemma 4.1. The edge grafting operation) Let u be a vertex of a graph G. Denote G u (a,b) the graph 

obtained by attaching to G two (new) pendent paths of lengths a and b at u. Let a,b,c,d be nonnegative 
integers with a + b = c + d. Assume that < a < b, < c < d and a < c. If u is a non-isolated vertex of a 
unicyclic or bipartite graph G, then the following statements are true: 

(1) . If a is even, then G u (a,b) >~ G u (c,d). 

(2) . If a is odd, then G u (a,b) -< G u (c,d). 

If a = 0, then we say that G u (0, b) is obtained from G u (c, d) by a total edge grafting operation. 
The following result in Q was obtained directly by using the edge grafting operation. 

8 



Theorem 4.1. fM] Let T be a tree of order n with N 3 (T) > 2. Then there exists a tree T' of order n with 
N 3 (T") = N 3 (T) - 1 and A(T') = A(T) suc/i f/iaf T -<! T'. 

In the followings, we will give some upper bounds for the energies of the trees of the form T n (a, b\c, d). First 
we consider the case 1 £ {a, b, c, d} in the following Theorem 14.21 The other case where min{a,b,c,d} > 2 
will be considered in Lemma 14.31 14.41 and Theorem 14.31 

Theorem 4.2. Q/ Let T = T„(l, b\c, d). Then T ~i P n (l, 2,n — 4). 

Proof. By using total edge grafting on the two pendent paths of lengths c and d, we have T -< P„(l, b, n—2 — b). 
Using the edge grafting operation again, we have P n (l, b, n — 2 — b) ^ P n (l,2,n — 4). Thus the result 
follows. □ 



Pa 




Fig. 6: G u , v (a,b) 



The following Lemma generalizes Lemma [4~T1 and is called "edge grafting operation at different vertices". 

Lemma 4.2. \1 w Let u,v be two vertices of a unicyclic or bipartite graph G. Let G UjV (a,b) be the graph 
obtained from G by attaching a pendent path of length a to u and attaching a pendent path of length b to v 
(as shown in Fig\$j). Suppose that G satisfies: 

(1,1)- 

(ii). For any nonnegative integers p,q, G u ^ v (p,q) = G u ^ v {q,p). 
Let a, b, c, d be nonnegative integers with a < b, c < d, a + b = c + d, and a < c, then we have 

(1) If a is even, then G u . v (a,b) >- G UiV (c,d). 

(2) If a is odd, then G u<v (a,b) -< G UtV (c,d). 

Fig. 7: Ti 2 (3, 2|2, 2) and T 12 (2, 2|2, 2) 

Now we use the methods given in §3 to prove the following two lemmas, which consider the tree T n (a, 2|2, 2) 
in two cases 3 < a < n — 9 and a = n — 8. These two lemmas will only be used in the proof of the Theorem 
Ol later. 

Lemma 4.3. Let3<a<n~9. Then T n (a, 2|2, 2) -< T„(2, 2|2, 2). 

Proof. Let ei,e 2 be the cut edges of G = Ti 2 (3,2|2,2) and e[, e' 2 be the cut edges of H = Ti 2 (2, 2|2, 2) 
as shown in FigJT] respectively. Then we have T n (a, 2|2, 2) = G(a — 3, n — 9 — a) and T n (2, 2|2, 2) = 
H(a — 3, n — 9 — a). 
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By some directly calculations, we have 

4>{H(0,0),x) =<£(r 12 (2,2|2,2),x) = x 12 + 11 x w + 43 a; 8 + 74 x 6 + 59 x 4 + 19 x 2 + 1, 
^(G(0,0),x) =<£(T 12 (3,2|2,2),x) = x 12 + 11 a; 10 + 43 a; 8 + 74 x 6 + 57 a; 4 + 17a; 2 , 
4>(H(l,Q),x) = 0(iJ(O,l),x) =0(T 13 (2,2|2,2),x) = a; 13 + 12 a; 11 +53 a; 9 + 108 x 7 + 107 x 5 + 48 a; 3 + 7 x, 
<£(G(l,0),x) =^(T 13 (4 ) 2|2,2),x) = a; 13 + 12 a; 11 + 53 a; 9 + 108 x 7 + 105 x 5 + 46 a; 3 + 7x, 
?(G(0,l),x) =0(7i 3 (3,2|2,2),x) = a; 13 + 12 a; 11 + 53 a; 9 + 108 x 7 + 106 a; 5 + 46 a; 3 + 6 x, 
0(77(1,1), a:) =0(Ti 4 (2,2|2,2),a;) = a; 14 + 13 a; 12 + 64 a; 10 + 151 a; 8 + 181 a; 6 + 107 a; 4 + 26 a; 2 + 1, 
4>(G(l,l),x) =0(7i 4 (4,2|2,2),x) = a; 14 + 13 a; 12 + 64 a; 10 + 151 a; 8 + 180 x 6 + 105 a; 4 + 25 a; 2 + 1. 

By comparing the coefficients of above polynomials, we find that 

G(0,0) -< if (0,0), G(0, 1) -< 77 (0, 1), G(1,0H 77(1,0), G(l, 1) -< 77(1, 1). 
So by TheoremElwe have T„(a,2|2,2) = G(a - 3, n - 9 - a) -<! 77(a-3,n- 9 - a) = T„(2, 2|2, 2). □ 

Fig. 8: G = Tu(3,2|2,2) and H = T u (2,2|2,2) 

Now we consider the remaining case a = n — 8 for the trees of the form T n (a, 2|2, 2). 
Lemma 4.4. E(T n (n - 8,2|2,2)) < E(T n (2, 2|2, 2)) /or a// n > 11. 

Froo/. Consider the cut edges e of G = T n (3, 2(2, 2) and e' of 77 = T xl (2, 2(2, 2) as shown in FigU Let G(fc), 
77(/c) be graphs obtained by subdividing the cut edges e of G and e' of 77 respectively k times. Then we have 
T„(n-8,2|2,2) = G(n - 11) and 7^(2,2,2,2) = H(n- 11). Denote 5fc = 4>{G{k),x) and 7i fc = 4>{H{k) lX ). 
By some directly calculations, we have 

fro =^(Tn(2,2|2,2),a;) = a; 11 + 10 x 9 + 34 x 7 + 48 a; 5 + 29 a; 3 + 6 x, 

9a =0(T u (3,2|2,2),a;) = a; 11 + 10 x 9 + 34 x 7 + 49 x 5 + 29 a; 3 + 5 x, 

hi =0(T 12 (2, 2|2, 2), x) = a; 12 + 11 x 10 + 43 x 8 + 74 x 6 + 59 x 4 + 19 x 2 + 1, 

9l =^(T 12 (4, 2(2, 2), x) = x 12 + 11 x 10 + 43 x 8 + 75 x 6 + 59 x 4 + 18 x 2 + 1. 

So we have 

Mo - h gi = x(x - l)(x + l)(x 6 + 7x 4 + 11 x 2 + l)(x 2 + l) 3 . 

Thus 

D = {x\h igo - h o9l < 0, x > 0} = (0, 1). 
Also by using computer we can find: 

E(77(0)) = 13.059967, E(G(0)) = 13.015698 
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and by using computer to calculate the integral we can further obtain 

1 

E(H(0)) - E(G(0)) + - /in ^rldx = E(H) - E(G) + - [ In ^dx = 0.005951 > 0. 
n J d (x) 7T J h gi 

D 

So using Theorem [231 we obtain E(H(k)) - E(G(k)) > for all k > 0. Thus E(T„(ra - 8, 2|2, 2)) < 
E(T n (2,2|2,2)). □ 

Theorem 4.3. Let n > 11, a, b, c, d > 2 and a, 6, c, d are noi a/Z egita^ to 2. Then we have 

E(T n (a,6|c,d))<E(T„(2,2|2,2)). 

Proof. By using the edge grafting operation in Lemma 14. 1[ we have 

T n (a,b\c,d) =4 T„(a + 6-2,2|2,c + d-2). 

By using Lemma T4.2I fedge grafting on different vertices), we also have 

T n (a + 6-2,2|2,c + d-2) =^ T n (a + b + c + d - 6, 2|2, 2). 

Write x = a + b + c + d — 6, then we have 3 < x < n — 8 since at least one of a, 6, c, d is greater than 2. 

Now If 3 < x < n - 9, then by Lemma [O] we have T„(a;,2|2,2) -<! T n (2,2|2,2). So E(T n (a,b\c, d)) < 
E(T n (x,2\2,2)) <E(T„(2,2|2,2)). 

If x = n-8, then by LemmaO] we have E(T„(a, 6|c, d)) < E(T n (x, 2|2, 2)) < E(T n (2, 2|2, 2)). □ 

5. The trees of order n with the first |_ "~ 7 ] largest energies 

In this section, we will determine the first L^r^J largest energy trees of order n > 31 by using the method 
of directly comparing energies given in §3. 

First, we divide the class of star like trees into the following four subclasses: 

(CI). The path P n . 

(C2). The class S n = {P„(2, a, b) \ a + b = n - 3, 1 < a < b}. 

(C3). The starlike trees T of order n with A(T) = 3 and T £ S n . 

(C4). The starlike trees T of order n with A(T) > 4. 

For convenience, we also define the following class (C5): 

(C5). The class of non-starlike trees of order n (i.e., N^(T) > 2). 

It is obvious that the union of the classes (C1)-(C5) is the class of all the trees of order n. 

Now, our strategy of proving the main result is as follows. Firstly, using the quasi-order we can obtain 
(in Theorem 1 5. ip a total ordering of all the L^y^J trees in S n . Secondly, we can show (in Theorem I5.2|) that 
the maximal tree (under the quasi-order) in the class (C3) is P„(4,4, *), and the maximal tree in the class 
(C4) is P n (2, 2, 2, *). Next, by directly comparing the energies of the largest energy trees in the classes (C3) 
and (C4) with some smaller energy graphs in S n , and comparing the energies of the tree T n (2, 2|2, 2) in the 
class (C5) with the smallest energy tree P n (2, l,n — 4) in S n , we obtain that the first L^Y^J largest energy 
trees in S„ together with P n are the first L^y^J largest energy trees in the class of all trees of order n. 
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Theorem 5.1. Let S n = {P n (2,a,b) \ a + b = n - 3, 1 < a < b}. Let k = L^J, t = [|J and I = [^-\. 
Then we have the following totally quasi order for the trees in S n : 

P„(2,2,*) y P n (2,4,*) y ■■■ y P„(2,2*,*) y P„(2,2Z + 1,*) >-•.•>- P„(2,3,*) >- P„(2,l,*). (5.1) 

Proof. The result follows directly from Lemma l4.il by using the edge grafting operation. □ 

Theorem 5.2. Let n > 11. T/ien we /lave 

(%). IfTe(C3) andT ^ P n (4,4,n-9), iAen T -< P n (4, 4, n - 9) . 

IfT e(C'4) andT ^ P n (2, 2, 2, n - 7), tfien T -< P„(2, 2, 2, n - 7) . 

Proof. (1) Since T €(C3), T must be of the form P n (a, 6, c) with 2 ^ {a, b, c}. Without loss of generality, 
we may assume that a < b < c. Then b + c > 7 since n > 11. So by Lemma 14.11 we have T = P n (a, b, c) =^ 
P„(a,4, 6 + c-4) and P n (a,4,6+c-4) ^ P„(4,4, n- 9) since fr + c-4 7^ 2. Also T 7^ P n (4,4,n-9) implies 
at least one of the above two relations is strict. Thus we have T = P n (a, b, c) -< P„(4, 4, n — 9). 

(2) Since A(T) > 4 for T S(C4), by using Lemma 14. II we can derive that T =4 Pn(a, b, c, d) for some tree 
P n (a,b,c,d). By further using the edge grafting operations at most 3 times on P n (a,b,c,d), we will finally 
obtain P n (a, b, c, d) =4 P«(2, 2, 2, n — 7). Also T 7^ P n (2, 2, 2, n — 7) implies at least one of the above relations 
is strict. Thus we have T -< P n (2, 2, 2, n - 7). □ 

The following Theorem 15.31 and Theorem 15.41 will exclude out P n (2,2, 2, *) (the maximal energy tree in 
the class (C4)) and T n (2, 2|2, 2) (in (C5)) by the smallest energy tree in S n by using the method of directly 
comparing energies given in §3. 




Fig. 9: P 9 (2, 2, 2, 2) and P 9 (2, 1, 5) 



Theorem 5.3. Let n > 10. Then we have E(P„(2, 2, 2, ra - 7)) < E(P n (2, 1, n - 4)) 

Proo/. Consider the cut edges e of G = P 9 (2, 2, 2, 2) and e' of P = P 9 (2, 1, 5) as shown in FigEl 

Let G(fc), P(fc) be graphs obtained by subdividing the cut edges e of G and e' of P respectively k times. 

Then we have P n (2, 2, 2, n - 7) = G(n - 9) and P n (2, l,n-4) = H{n-9). Denote = fi(G(k),x) and 

h k =$(H(k),x). 

By some directly calculations, we have 

/i =0(P 9 (2, 1, 5), x) = x 9 + 8x 7 + 20x 5 + 17x 3 + 4x, 
go =0(P 9 (2, 2, 2, 2), x) = x 9 + 8 x 7 + 18 x 5 + 16 x 3 + 5 x, 

hi =0(Pio(2,l,6),x) = x 10 +9x 8 + 27x 6 + 31x 4 + 12x 2 + l, 
3i =0(Pio(2,2,2,3),x) = x 10 + 9x 8 + 25x 6 + 28x 4 + 12x 2 + l. 

So we have hig - h gi = (2 x 4 + 8 x 2 + l)(x 2 + l) 3 > for all x > 0. 

Also we can compute that E(P(0)) = E(G(0)) = 6 + 2 V5. So using Theorem 13. II we have 

E(P„(2,l,n-4)) -E(P„(2,2,2,n- 7)) = E(H(n - 9)) - E(G(n - 9)) > E(P(0)) - E(G(0)) =0. □ 
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Notice that P n (2, 2, 2, n — 7) and P„(2, 1, n — 4) are quasi-order incomparable when n > 11. So Theorem 
15.31 can not be proven by only using the quasi-order method. 




Fig. 10: T 22 (2,2|2,2) and P 22 (2, 1,18) 



Theorem 5.4. Let n > 22. Tften we have E(T n (2, 2|2, 2)) < E(P n (2, 1, n - 4)). 

Proof. Consider the cut edges e of G = T 22 (2, 2|2, 2) and e' of H = P 22 (2, 1, 18) as shown in FigOU 

Let G(k), H{k) be graphs obtained by subdividing the cut edges e of G and e' of H respectively k times. 

Then we have T n (2,2|2,2) = G(n - 22) and P n (2,l,n - 4) = H(n ~ 22). Denote g k = 4>{G{k),x) and 

h k = j(H(k),x). 

By some directly calculations, we have 

h = x 22 + 21 x 20 + 189 x 18 + 953 x w + 2955 x u + 5824 x 12 + 7293 a; 10 + 5643 x s + 2541 x 6 + 595 x 4 + 57 x 2 + 1, 
3o = x 22 + 21 a; 20 + 188 x 18 + 939 x 16 + 2879 x 14 + 5625 x 12 + 7046 x 10 + 5546 x s + 2598 x 6 + 644 x A + 64 x 2 + 1, 
ftx = x 23 + 22 x 21 + 209 x 19 + 1123 x 17 + 37 56 x 15 + 8113 x 13 + 113 75 x 11 + 10 1 53 x 9 + 5511 x 7 + 16 72 x 5 
+ 241x 3 + 11 x, 

,9i = x 23 + 22 x 21 + 208 x 19 + 1108 x 17 + 36 67 x 15 + 78 50 x 13 + 109 82 x 11 + 99 1 2 x 9 + 55 46 x 7 + 1768 x 5 
+ 268x 3 + 12 x. 

So we have 

h l9o - h Q gx = x(x 8 + 7x 6 + llx 4 -4x 2 - l)(x 2 + l) 3 

D = {x\hxg Q - h g! < 0,x > 0} = (0,0.663073). 
By using computer we can also find 

E(H(0)) = 27.182092, E(G(0)) = 27.175139, and E(P(0)) - E(G(0)) + - f ln(^)dx = 0.000425 > 0. 

7T J h gx 

D 

So by using Theorem [XI we have E(P„(2, 1, n-4)) -E(T„(2, 2|2, 2)) = E(H(n- 22)) -E(G(n- 22) > 0. □ 




Fig. 11: P 3 i(4,4, 22) and P 31 (2, 7, 21) 



The following Theorem 15.51 will exclude out the maximal energy tree in the class (C3) by the fourth 
smallest energy tree in S n - 
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Theorem 5.5. Let n> 31. T/ien we have E(P n (4, 4, n - 9)) < E(P„(2, 7, n - 10)). 

Proof. Consider the cut edges e of G = P 3 i(4, 4, 22) and e' of PT = P 3 i(2, 7, 21) as shown in FiglTTl 

Let G(fc), H{k) be graphs obtained by subdividing the cut edges e of G and e' of P respectively times. 

Then we have P„(4, 4, n - 9) = G(n - 31) and P n (2, 7, n - 10) = iJ(n - 31). Denote o fe = 0(G(fc), a;) and 

h k = 0(H(k),x). 

By some directly calculations, we have 

ho = 0(P 3 i(2,7,21),x) = x 31 + 30 a; 29 + 405 a: 27 + 3252 x 25 + 17296 a; 23 + 64 2 20 x 21 + 170 9 43 a; 19 + 329768 x 17 

+460696 a; 15 + 460851 x 13 + 322620 a; 11 + 152131 a; 9 + 45426 x 7 + 7738 x 5 + 619 x 3 + 15 x, 
go = <?KP 3 i(4,4,22),x) = a; 31 + 30 a; 29 + 405 a; 27 + 3252 a; 25 + 172 95 a; 23 + 64 2 00 a; 21 + 170772 a; 19 + 328952 a; 17 

+458317 a; 15 + 456496 a; 13 + 317681 x 11 + 148864 a; 9 + 44349 a; 7 + 7644 x 5 + 636 a; 3 + 16 x, 
^ = 0(P 32 (2, 7, 22), x) = x 32 + 31 a; 30 + 434 a; 28 + 36 29 a; 26 + 20 1 98 a; 24 + 78938 a; 22 + 22 2 7 24 a; 20 + 45 9 3 65 a; 18 

+693530 x 16 +760145 a; 14 + 593801 a; 12 + 320464 a; 10 + 113705 x 8 + 24470 x 6 + 2774 x 4 + 125 x 2 + 1, 
9i = 0(P 32 (4, 4, 23), x) = x 32 + 31 a; 30 + 434 a; 28 + 36 29 a; 26 + 20 1 97 a; 24 + 789 1 7 x 22 + 222534 a; 20 + 45 8 3 96 a; 18 

+690471 x 16 +753971 x 14 + 585871 x 12 + 314249 x 10 + 111032 x 8 + 24007 x 6 + 2792 x 4 + 132 x 2 + 1. 

So we have 

Mo ~ h g! = x(x 4 + 3x 2 + l) (x 12 + 12 x 10 + 53x 8 + 107 x 6 + 99x 4 + 34x 2 + l) > for all x > 0. 
By using computer we can also find 

E(iT(0)) = 38.616923, E(G(0)) = 38.616742 

So using Theorem O we have E(P n (2,7,n- 10)) - E(P„(4,4, n - 9)) = E{H{n - 31)) - E(G(n - 31)) > 
E(H(0)) - E(G(0)) = 0.000181 > 0. □ 

Theorem 5.6. Let n > 31. Let S' n = S n \{P n (2, 5, n - 8), P„(2, 3, n - 6), P„(2, 1, n - 4)} be the first [^\ 
trees in the quasi-order list i5.1)) of S n . Then P„ and the [-^-^J trees in S' n are the first L :2 x^J largest energy 
trees in the class of all trees of order n. 

Proof. It is obvious by the quasi-order list (|5.1I) that the smallest energy tree in the set {P n }L)S' n is P n (2, 7, n— 

10). Now take any tree T ^ {P n } U S' n of order n, we consider the following four cases: 

Case 1: T e(C2). Then T e S n \S' n . By the quasi-order list we have T -< P„(2, 7, n - 10). 

Case 2: T e(C3). Then by Theorem [O] and Theorem [53] we have 

E(T) < E(P„(4, 4, n — 9)) < E(P„(2, 7,n— 10)). 

Case 3: T G(C4). Then by Theorem O E31 and the list (JSTTJ) we have 

E(T) <E(P„(2,2,2,n-7)) < E(P„(2, 1, n - 4)) < E(P„(2, 7, n - 10)). 

Case 4: T e(C5). 

Subcase 4.1: N 3 (T) = 2 and A(T) = 3. Then T is of the form T n (a,b\c,d). So by Theorem E21 |43J E3] 
and the list (|5.ip we have 

E(T) < E(P„(2, 1, n - 4)) < E(P„(2, 7,n - 10)). 
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Subcase 4.2: N 3 (T) = 2 and A(T) > 4. Then a tree T with i\T 3 (T') = 2 and A(T') = 3 can be obtained 
from T by using total edge grafting several times. So T -< T", and thus by Subcase 4.1 we have E(T) < 
E(T') <E(P„(2,7,»-10)). 

Subcase 4.3: Ns(T) > 3. Using Theorem 14.11 several times we can obtain a tree T' with A^(T') = 2 and 
T -< T. So by Subcases 4.1 and 4.2 we have E(T) < E(T') < E(P n (2, 7,rc- 10)). □ 
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